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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

A SPECIAL METHOD FOR FINDING BODY DISTORTIONS THAT 
REDUCE THE WAVE DRAG OF WING AND BODY 
COMBINATIONS AT SUPERSONIC SPEEDS 
By Harvard Lomax and Max. A. Heaslet 

SUMMARY 


For a given wing and supersonic Mach number^ the problem^ of shaping 
an adjoining fuselage so that the combination will have a low wave drag 
is considered. Only fuselages that can be simulated by singularities 
(multipoles) distributed along the body axis are studied. However^ the 
optimum variations of such singularities are completely specified in 
terms of the given wing geometry. An application is made to an elliptic 
wing having a biconvex section^ a thickness -chord ratio equal to O.O 5 at 
the root^ and an aspect ratio equal to 3- A comparison of the theoretical 
results with a wind-tunnel experiment is also presented. 


INTRODUCTION 


The most simplifying assumptions that still permit the construction 
of a mathematical model general enou^ to contain quantitative informa- 
tion about steady three-dimensional supersonic flow are those basic to 
the development of linearized theory. Of these ^ the two principal assump 
tions are that the viscosity effects are negligible and the perturbation 
velocities are almost everywhere small enou^ to be neglected relative 
to the fli^t or free-stream velocity. Under such restrictions the flow 
field can be described in terms of a perturbation velocity potential cp 
obeying the equation 




( 1 ) 
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where 


1 and the reference coordinate system^ is shown in 

sketch (a). Further, the wave drag 
of any object in a flow field governed 
by equation (l) can be evaluated (see, 
e.g. , ref. l) by means of the equation 



D = 


-“of 


2 Jt r 
d0^ 


lim (rcp^cp^) | (2) 


00 


'] 


where x,r, and 0 are cylindrical 
coordinates also defined in sketch (a). 

General solutions to equation (l). 

. are numerous and classical. In apply- 
ing these solutions to the interpreta- 
tion of physical phenomena the usual 
approach is to fit them to the given 
boundary conditions, that is, to make 
the flow field simulated by them con- 
form to the shape of the disturbing 
object as well as to a uniform free 

stream at infinity. Hence, from this point of view, the choice of a type 
of general solution to be used in analyzing a particular problem with the 
least mathematical effort depends on the geometric form of the object 
under consideration. For example, general solutions based on Green's 
theorem are well adapted to the study of forces on single planar wings 
in a steady supersonic flow. On the other hand, the general solution^ 
given by Lamb (ref. 2) - which is composed of an infinite set of multi- 
pole distributions disposed along a line - is well adapted to the study 
of the flow around fuselage-like objects. 


In this report use is made of certain general solutions to equa- 
tion (1) but with a deviation from the usual approach mentioned above. 

One considers, in fact, two different kinds of solutions which represent 
separately, in a given vicinity, different classes of real objects ^d, 
by means of equation (2), finds optimum combinations of these solutions 
from the viewpoint of low wave drag. The analysis involved in solving 
this problem has, in general, a distinct mathematical advantage over the 
problem of calculating the drag of a given object; namely, that the 
immediate problem of finding a shape with a relatively low wave drag is 
divorced from any detailed reference to the shape itself. It is true, 
of course, that the stream surface representing this shape must eventually 
be found and, in fact, a limitation on the applicability of the method 
is given by the requirement that this shape be real. However, the problem 
ijt should be stressed that the x axis is parallel to the free- 
stream direction (wind axes) so a body of revolution can be symmetrical 
about this axis only at zero angle of attack. 
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of finding the shape of the object when cp is known is a matter of direct 
calculation.^ One should also be careful to notice that the optimum solu- 
tions obtained by this procedure are not necessarily true optimums but 
purely relative to the choice of solutions used in the analysis. 


LIST OF IMPORTAOT SYMBOLS 


A 

An(x) 


a 

ajj( x) 


Bn(x) 


b 

^n 

Cd 

Cp 


wing aspect ratio 

times the nth derivative of the nth multipole dis- 
tribution an(x) (See eq. l6.) 

semiroot chord of elliptic wing 

strength of nth-order multipole distribution multiplying 
cos ne 

(-p)’^ times the nth derivative of the nth multipole dis- 
tribution bn(x) 

semispan of elliptic wing 

strength of nth-order multipole distribution multiplying 
sin ne 

drag coefficient, 

pressure coefficient, local pressure minus static pressure 
divided by q 


D wave drag 

Dji wave drag associated with nth-order multipole distribution 

(See eq. 59*) 

Lo*,Lq maximum fore-and-aft extent of wing equivalent multipole 

distribution 


l’( 0),L(0) maximum fore-and-aft extent of wing equivalent multipole 
distribution for angle 0 


M free -stream Mach number 

^From a mathematical point of view the essence of the method outlined 
above is that the analysis involves the solution to direct problems, that 
is, problems of integration. Calculating the drag of a given body, on the 
other hand, involves the solution to inverse problems, that is, problems 
involving the inversion of integral equations^ 
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q 

f 


fe 

R 

S 

S^(x,6) 

t 


Uo 

V 


x,y,z 

X,r,e 


a(x) 


, P 

Po 

9 


PoUo^ 

free-stream dynamic pressure, — g — 

E 

Lo 

see equation (46) 

radius of tody 

area of wing plan form 

normal projection of wing cross-sectional area measured in 
oblique planes 

maximum thickness of wing root chord 

speed of free stream 

volume 

Cartesian coordinate system, x parallel to free-stream 
direction 

cylindrical coordinate system, x parallel to free-stream 
direction 

strength of nth-order optimum cancellation multipoles 
(See eq. (36) .) 

^/M^ - 1 

slope of wing upper surface measured parallel to free-stream 
direction 

tan ^ (P cos 0) 

free-stream density 

perturbation velocity potential 

DEFINITION OF THE PROBLEM 


The problem of designing an airplane to have a minimum wave drag 
must be stated quite precisely. If the aerodynamic is t is approached with 
the question, "Given an aerodynamic shape, can its wave drag be lowered? 
he can always reply that any volume of material having a wave drag can 
always be reshaped within a space of finite dimensions so that it will have 
less wave drag at a given Mach number. Such an answer is interesting but, 
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at present, not very useful to the airplane designer. There is first, 
of course, the basic criterion that the total drag should be minimized 
at a given lift and minimizing a component part of this total without 
holding the other parts fixed does not necessarily yield the lowest pos- 
sible drag for a given set of restraints. For example, the configuration 
illustrated in sketch (b) has no wave drag when traveling at zero angle 
of attack; but it has a relatively 


Cylindrical shroud 


hi^ friction drag, because of the 
large amount of wetted area, and its 
drag due to lift could also be rel- 
atively high. Completely aside from 
all such performance considerations , 
however, are many other important con- 
siderations that are unfortunately more 
or less vaguely defined from an aero- 
dynamic point of view. For example, 
an airplane must contain a certain 
amount of usable volume, the shaping 
of individual parts is limited by 
structural requirements, and the 
arrangement of these parts must not 
seriously harm the airplane stability 
and control. The interrelation of 
all such separate demands presents an 
extremely complex design problem making it difficult to deviate too far 
from the reliable shapes set by experience. 



— Mach waves 


Sketch (b) 


As a result of the above-mentioned difficulties, the aerodynamicist 
who is concerned with discovering a practical airplane shape having low 
wave drag finds the real definition of his problem somewhat obscure. In 
a sense his first problem is, literally, to pose a problem; that is, to 
impose a minimum number of arbitrary but pertinent restraints within the 
framework of which the wave drag is to be minimized. Even when this has 
been done, he still is concerned with the question of uniqueness, since 
optimum shapes are not necessarily unique even when several restraints 
are imposed. Consider, for example, the problem of finding the Busemann 
biplane which will have minim\mi wave drag at a given Mach number for a 
fixed section strength, volume, and 
wetted area. If the design Mach 
number is 1.4l, one such design (on 
the basis of linearized theory) is 
shown in sketch ( c) where the chord- 
gap ratio, h/c, is equal to l/2. The 
resulting variation of the wave drag 
is shown in the upper part of sketch 
(d). However, when the gap is closed 
to the point where h/c equals l/4, 
the variAtion of wave drag, shown in 


Mach waves at M = 1.41 

\ ’a 




f- 


c- 

Sketch (c) 


/T 

h 
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/ 1.2 1.4 1.6 t.8 2D 

M 



/ L2 1.4 t.6 1.8 2D 


M 


Sketch (d) 


the lower part of sketch (d), is 
the same within the interval 
1.28 < M < 1.66 and everywhere 
else is lower. It is likely that 
one would have first discovered 
the former solution, yet to the 
accuracy of the theory used, the 
latter is obviously preferable. 

With the above observations 
always in mind, attention will be 
directed in this -report to the 
analysis of simplified configura- 
tions composed of two distinct 
types of volume: planar types, 

that is, wing-like volumes, thin 
in one dimension and bounded by 
surfaces that never deviate far 
from a reference plane; and rec- 
tilinear types, that is, fuselage- 
like volumes longer in one dimen- 
sion than in the other two and 
disposed more or less symmetrically 
about a straight line 

In particular the following 
problem is posed: 


Given a thin nonlifting wing, what is the shape of 
an adjoining fuselage, the stream surface of which 
is simulated by a line of multipoles in the same 
plane as the wing, that will minimize the wave drag 
of the combination at a given Mach number? 


BASIC CONCEPTS 
A Line of Sources 


The velocity potential induced at the point x,r,0 by a group of 
sources distributed along the x axis, starting at -Lq, is well known 
to be given by 

ao(|)d| 




cp(x,r,e) = - 


n/(x - 1)2 - p2j.2 


(3) 
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where ao(|) is the source strength per unit of length. In order to 
calculate wave drag one needs only the value of cp as r approaches 
infinity. This asymptotic value is simple enough to find provided it 
is observed that, as r ' is increased, x should also he increased so 
the potential can he studied in the vicinity of the Mach waves radiating 
from the disturbing object. Hence, set 


X = Xo + Pr 


W 


so for a given r, Xq measures the streamwise distance of the point 
x,r,e from the Mach wave emanating from the origin and,, in particular, 
the foremost wave is located at Xq = "Lq. (See sketch (e).) 

If equation (4) is placed into 
equation ( 3 ) and r is assumed to 
he large, the potential induced by 
a source is 




<p) 


1 ap(|)d| 


(5) 


and the induced velo cites are 


fPx) 


'x^r— ^ 00 


; /’ ° ao*(Qd| 

J - n/xo - | 

“Ln 


(6a) 


'x^/3r / 


P(x,r,$) 


\i \'' 




■/Sr 


Sketch (e) 


Xf,-*\ 


(pj.) 


r'r-^ 00 


= -Pcpi 


(6h) 

Multipoles 


Lamb, in reference 2, page 527, has presented a general solution to 
equation (l) consisting of an infinite set of basic singular solutions. 
These basic singularities, referred to as multipoles, can be distributed 
along a line and wei^ted so as to reproduce certain body shapes enclosing 
the line. The expression for the perturbation velocity potential for a 
distribution of nth-order (n =0,1,2, ...) multipoles starting at -Lq 
and continuing along th® x axis can be written in terms of a cylindrical 
coordinate system (sketch (a)) as 
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/ A S '? [anU)cos n0 + bn(|)sin ne]d| ^ ^ 

^n(x,r,e) - - 5, (, 5j; / 


The operator ^ ) is defined as 


1 sY =1 s 4 = - i ^ 

r Sry r ^ \r r® ^ r® Sr^ 


and the definition of (- ) follows by induction. If the notation® 

V* Sr/ 


A(y)dy , (1)(3) ■ — (2n - 1) 


= (- 1 ) 


'X - y 


A(y)dy 


(X - y) 


( 2 n+i)/s 


is introduced, where the symbol f is read "finite part of the integral," 
equation (7) becomes 


r“p®“(2n)! [an(|)cos ne + bn(l)sin nO] d| 

cpn(x,r,0.) = - . n+i , , , 7 (pn-fi)/ 7~ 


and the general expressions for the induced velocities become - writing 
only the term involving the cosine, since the result for the sine is 
identical - 


cpx(x 


oo 

> V r°H^°(2n g)! r 
’ ’ L 2="+=„(n + 1)! 7 


(x - |)an(|)cos n0 d0 


Lo [(x - 1)^ - P^r^] 


2 o 2 2l(2n+3)/2 


V P^’^r"“^(2n) ! P^'^^ [n(x - (n+ l)P®r^]an(5)cos n0 d^ 

q> 3 ,(x,r, 0 ) = pn+i„(n) A .,2 . a,(2n+3)/2 


[(x-l) -pr ] 


(an+s) /2 


®For a detailed discussion of the finite-part concept as used in this 
report see reference 3* 
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1 , , V f 

■ Z. 2-^(n)T"‘J 


x-pr 


and) sin a© d| 


[(x - if - p2r2] 


( 2 n+i )/2 


(9c) 


Another very useful way of developing these multipole solutions 
evolves from an application of operational techniques. To begin with, 
rewrite equation (l) in terms of a polar coordinate system, thus 


2 ^ - 0 


^ ^ ^ 3r r^ 


( 10 ) 


Next, define the Laplace transform of cp(x,r,e) by 


^(s,r,e) = J' cp(x,r,0)e"^Sdx 


( 11 ) 


and apply this transform to equation (lO). There results (for a proof 
see Appendix A) 


P^s‘ 




■cp - 


Sr 


2 _ 

± n. - — £_£ = 0 


Sep 

Sr 




( 12 ) 


Now, if a general solution to equation (12) is expressed in the form 

(p(s,r,0) = f(r)cos n0 
then f(r) must satisfy the equation 


+ 1 M _ + fs^)f=0 

dr^ r dr Vr2 J 


Solutions to this are given by 


f(r) = Cn(s)ln(Prs) + An(s)Kn(prs) 


where In ’and Kn are modified Bessel functions as defined in reference 4, 
page 77. Hence, if ^ is to vanish as r goes to infinity, a general 
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solution to equation (lO) can te written in the form 

00 

9(s,r,e) =^^[In(s)cos n© + Bn(s)sin n©]Kn(prs) (I3) 


The above result will he transformed hack to the physical plane in 
two ways. First, apply the identity (ref. k, p. 79) 

Kn(.) = (-1 )“(|^)ko(^) 

and re-express equation (I3) as (only the coefficient of the cos n0 
term is written since the treatment of the sine term is identical) 

00 ^ 

$(s,r,0) = - ^^(-l)^cos n0 r’^ Ko(prs) 


The inverse Laplace transform of Kg(prs) - see reference 5 " is 


L [Ko(prs)] = < 


0 

1 


i Vx2 - 


P^r® 


X < Pr 
X > Pr 


So, since 




^n-i 

dxnAn( xn) 


-L 


o 



n-i 

Xi) An(xi)dxi 


(11^) 


an application of the convolution integral and other standard operational 
techniques yield 


cp(x,r,0) = - 


L 


x-pr 


Ao(0<i| 


n/(x - 1)^ - p^r^ 


00 -I 

i-a&f 


n pX-pr 


(- 1 ) 


n 


n-i 


[?wT)T.L(' 


''/(x “ I )^ - P^r® 


(15) 
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From compgirison of equations (8) and (15) > the relation between the 
strengths an(x) and An(x) for the two different forms of the solution 
is found to be 


ao(x) = Ao(x) 


(-p)“an(x) 



n-i 


xi) An(xi)dxi, 


» ( l6a) 

n > 0 


or 


= An(x) 


(l6b) 


where a^’^^(x) symbolizes the operation £41(2^) where use is made 


of the conditions 


(n)/ T \ T ) 

3n ( “Bq) = Sji { -Lq) - 


— ajj(-LQ) = 0 


(17) 


Another way to transfonu equation (I3) back to the physical plane 
is to do so directly. In this way one finds (from ref. 5) 


0 


h [Kn(prs)] - I coshj^n cosh“^^^^ 


'/x^ - p^r^ 


X < Pr 


X > Pr 


from which equation (I3) reduces immediately to 


cp(x 


00 


An(|)cosh 


n cosh 


n=o 


n/(x - 1)^ - p\^ 




d| 


(18) 


The pertxirbation velocities in the field represented by this potential 
are readily calculated. Thus 
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9 x(' 


“ pX-pr An’ ( 0 coshln cosh" 

:,r,e) = - i) cos ne / ===== ^ 

L - 5 )^ - pV 


(19a) 


00 

icp0(x,r,0) = sin nej 


And) cosh 


n cosh 


i/x - I 


7 (x - 1)^ - 


(I9h) 


and hy taking the derivative of equation (I3) with respect to r, one 
finds 

00 > 

= -i |psAo(s)Ki(Prs) + | ^PsAn( s)[ Kn-i(Prs) + Kn+i(Prs)]j- 

1 


which transforms to 

1 f px-pr (x-e)V(e)ig ^ p 

rJ(x-if -pV 


wxixuxi ox euicj j. 1110 

\-Lc 


■Lo 

x-Pi* Ar,'fE^^cosh 


An'(l){c 


(n-l)cosh 


tli)] 


+cosh 


-Lr 


n/(x-|)^-PV 


(n+1) cosh"^ 


(19c) 


If the relation between the functions an(x) and An(x) is given by 
equation (16), the velocities represented by equations ( 9 ) and (I9) are, 
of course, identical. 

In order to obtain limiting values induced by multipoles distributed 
along the x axis starting at -Lq, one returns to either equation ( 8 ) 
or (18) and calculates the leading term in a l/r expansion. As in the 
derivation of equation ( 5 ), it is necessary to observe that as r^ is 
increased, x should also be increased so cp is given in the vicinity 
of the foremost Mach cone created by the multipole distributions. Hence, 
using equation ( 4 ), one finds for equation (18) 
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ij^n cosh ^(1 + 


^ “ Xo An(Ocosh| 

^.--^oosnej 

o -Lo y2pr(xo - |)(^1 + 


Xq - ( 

3r 


m 


Xq - i 
2|3r 


which has the leading term as r goes to infinity 


cp) 


r ^00 




^ cos J^. 
o -L^ 


AnU)dg 

n/ Xq “ r 


(20) 


Similarly, the perturbation velocities reduce to 


9x) 


0° I / V 

1 V r ° An U)d| 

) cos ne / 


U -ijQ 


(21a) 


^00 *^*^x 


(21b) 


r ^9 


r 00 


^ynslnner°^ 

2n^2pr Lm J n/xo~“ 

o “Lo 


U)£l 

n/xq “ I 


(21c) 


In calculating the vave drag using equation (2) only the velocity 
components 00 > 00 necessary. Hence, from comparison 

of equations (21a) and (21h) with (6a) and (6b), it follows that at a 
given 6 a series of multipoles" induce the same momentum flux on an 
infinite cylindrical control surface^as a line of sources having a 
strength variation ao(x) equal to ^ cos nOkj^ii) * If one identifies a 
line of sources with a body of revolution, then it is apparent that, at 
a given 0, a dragwise equivalence has been established between a line 
of multipoles and a body of revolution. 
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Hayes’ Theorem and Its Application 


In the previous section a relation was found between multipole and 
source strengths which produce, at a fixed 0, equivalent momentum trans- 
port across a cylinder of infinite radius. By using a theorem due to 
Hayes (ref. 6) one can derive the strength relationship between any dis- 
tribution of singularities throu^out space and a line of sources which 
gives the same equivalence. 

The essence of Hayes* theorem is that, for a fixed 6, the velocities 
induced on a cylinder of infinite radius by singular solutions to equa- 
tion (1) (e.g., sources and doublets) are invariant to displacements of 
the singularities along certain oblique planes. In order to be specific, 
the equation of these oblique planes is next derived. 

Consider the point x,r,0 in a flow field having a supersonic free 
stream moving parllel to the x axis. Sketch (f) shows the Mach fore- 
cone (by definition the Mach forecone is the boundary of the region within 
which a disturbance in a supersonic stream can affect the flow at the 



cone apex) from x,r,0 in Xx^yi^Zi space. The equation of the forecone 

is 

xi = X -^^/(r cos 6 - yi)^ + (i* sin 0 - z^)^ (22) 
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One wishes to let r become very large and find the shape of the fore- 
cone as it passes throng regions close to the origin of the Xi,yx^Zi 

coordinate system, regions in which the objects creating the wave drag 
are located. From equation (4), and the expansion of equation (22) for 
large r, it follows that 

xi = Xq + pr - pr 1 - § (yicos 6 + zisin 0) + 

= Xq + P(yicos 0 + zisin 6) - ^ (yisin 0 - zicos 0)^ + ♦ • ♦ 

and when r goes to infinity 

Xi = Xo + Pyicos 0 + pz^sin 0 (23) 

which is the equation of the oblique plane mentioned above. It should 
be noted that the envelope formed from these planes by fixing xq and 
varying 0 between 0 and 2it coincides with the Mach forecone and after- 
cone from the point Xq,0,0. 

4 

Hayes’ result can now be stated 

To the lowest order. in l/r, as r tends to infinity, 
the magnitude of the perturbation velocity potential 
and its gradients at a fixed azimuth angle is invar- 
iant to a finite translation of sources (or any other 
singular solution to the wave equation) on planes 
parallel to that given by equation (23). 

Consider the application of Hayes' theorem to planar distributions 
of sources lying in the = 0 plane. As is well known, such a distri- 
bution simulates a wing symmetrically disposed about the horizontal 
(zi = 0) plane. In fact, if Au(xi,yi) is the local slope of the wing 
upper surface, the local source strength per unit area (according to thin 
airfoil theory) required to simulate the wing is -Uo^u/jf and the veloc- 
ity potential of the disturbed flow field is given by 



\i(xi,yi)dxidyi 

^(x - xi)^ - p^(y - yi)^ - P^z^ 


^For proofs, see Hayes' original derivation (ref. 6) or, if more 
convenient, reference 1. 
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where t is the area of integration 
(see sketch (g)) bounded by the wing 
edge and the trace in the = 0 
plane of the Mach forecone from the 
point x,y,z. Next introduce, the new 
coordinates |i and t]i such that 
lies along the Xi axis and t^x lies 
along the intersection of the Zx = 0 
plane and the plane given by equa- 
tion (23) (see sketch (h)). Set 


tan n = P cos 6 (25) 


and 

= xx - yxian 


■Hi = yisec (I 


V (26) 


xi = + Tixsin u 


yx = mcos |i J 


Then, in terms of the Ii,tix system, 
equation (2k) becomes 


cp(x,y,z) = - ^ r f — 

T Ji 


Mil + Tlisin jjL,TixCOs ti)cos [i dlxdTix 


X - li - ilisin p)‘ 


P^(y - mcos n)^ - P^z^ 


(27) 


As before, the asymptotic value of cp as r = -7 y^ + z^ » 00 is to be cal- 
culated. Accordingly, one can apply Hayes* theorem and sum up all the 
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sources along a line = constant (e.g.^ between a. and b in sketch (i)) 
and place them as a single source on the axis. The strength of this 
large equivalent source is given by 


UoSw - 2 UqCOS |i 


f ^u ( 1 1 


+ qisin |i,qiCOs ^l)dT^; 


(28) 



^The true oblique plane is given by equation ( 23 ) but the wing 


is so close to the = 0 plane that the variation with can be 


neglected. 
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no loading (local lift) will be considered. Lifting effects can be 
treated in a similar fashion, of course , if the wing doublet strengths 
are known. 

Cancellation Multipoles and Drag Minimization 

Since the flow field is governed by a linear partial differential 
equation the velocities induced by different solutions to it are additive. 
Therefore, the drag of an object simulated by various multipoles dis- 
tributed along the Xx axis and a sheet of sources in the Zi = 0 plane 
is given by ■ 


D = -Po 


r-iH 


lim 
r 00 


(9r)m + (^r) 


w 


(<Px)m 


w 


where the subscripts m and w refer to the multipoles and wing sources, 
respectively. But equations (21) and (29) identify, for a fixed 0, 
these velocities with those induced by equivalent line sources. Hence, 
for any given 0, oiie can immediately apply Karman*s drag formula (ref. 7) 
and then for the total drag, integrate 0 from 0 to 2it. This leads to 


D = - 


_Pp 

8rt' 


pSTi pL(0) ^L(0) 

^ / d0 / dxi / dx2 
' " -l’(0) -L'(e) 


00 

UoSw"(xi,e)+^(-B)“|^a^“'^^^(xi)cos n0 + bn (xi)sin n0 
" o 

00 . \ 
|uoSw”(x2,0)+^(-3)“|^a4“'^^^X2)cos n 0 + bn^"^^ (xg) 


)sin n0 


j-lnlxi- X2 


o 

However, since both the wing and multipoles are in the same plane (inter- 
preted physically, the wing is centrally mounted on the fuselage) and 
the wing is simulated by sources only (has no twist or camber), one can 

show the optimum value of each (xx) is identically zero.® 


6 By symmetry 


Hence 


Sw"(x,0) = Sv”(x,2rt - 0) 

/ 2J( 

sin n0Sw"(x,0)d0 = 0 


b (n-Hi) 

and any positive or negative variation of bn (xx) can only increase 

the dbag. 
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Hence, one can write 


p pl(6) p^(o) r ^ n Tti-i.t'* 1 

J <10 / <1x2 UoS^"(xi,e)+2^ (-3) an ^(xi)cosn9 

o -L'(e) -l’( 0) ^ o -* 


00 

lJoSw"(x2,e)+y (-P)’^an’^'^^^(x2)cos n0 


U 

O 


ln|xi - X 2 


(30) 


Next expand the term Sv"(x,0) in. a Fourier series. One finds 


UoSw"(x,e) 

O 


(x)cos n6 


(31) 


where 


u 

^o(^) ” 2 ^ j ^ (x^ 0 )d 0 


(32) 


<=<4^^(x) = — n / S^*(x,0)cos n0 d0 (33) 

n(-p) ^ 

Place these expressions in equation (3O), integrate with respect to 0 - 
using the orthogonal property of the trigonometric series - eind one finds 


D = 2 Dq + ^ 

1 


00 


where 


q 


2 n 


4 ltUn^ 


/ 

•"Ln "l»n 


dx- 


■ (n+l) (n+l) 

of-n (xx) + an (^i) 


(n+l). X (n+l) xl , I 

On (X2) + an (X2) ln|xi - X2 | 


( 34 ) 


( 35 ) 
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Since one can show”^ 


I = - / ±xi I dx2f’(xi)f»(x2)ln|xi - X 2 I >0 


-Ln “Lf 


the minimum value of D as expressed by equation (34) is given when 
each is itself a minimum. In other words, each Dn can be mini 

mized separately. Further, it follows that the value of the minimum 
itself is zero and occurs when 




Equation (36) is the mathematical definition of the optimum cancel- 
lation multipoles; namely, those multipoles which are just equal in mag- 
nitude and opposite in sign to the wing equivalent multipoles - equivalent 
in the sense that they induce an identical momentum flux across a cylinder 
of infinite radius. 

Obviously, if all the optimum cancellation multipoles were used, the 
wave drag of the combination would be zero. This result must, however, 
be properly interpreted with regard to the simulated shape. In order that 
the multipole lines can represent the distortion of a real fuselage, one 
must assume a cylindrical body exists upstream from the Mach cone 
X + Lo = pr (the effects of the nose are being neglected). This body 
forms the initial boundary of the stream tube which represents the physical 
fuselage in the vicinity of the wing and multipole lines. Clearly, the 
area enclosed by this initial boundary can be small enough for the sub-^ 
sequent stream surface to cross itself and represent, therefore, a physi- 
cally unreal body. Hence, the fact that the wave drag of the wing and 
multipole combin ation can be reduced to zero is quite valid, but in the 

’’’set f(x) = E^^sin ni, x = -Lqcos \lr. Integration gives 


00 



X 


. Ho 
2« J '■ 


Sw* (x,e)de, , n = 0 


2rt 

Ho__ r s "(x,e)cos ne de, n >0 

(36) 


which can never be negative. 
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over-all picture not only have the inevitable nose and tail drags been 
neglected but also the shape simulated by the combination can be unreal 
istic. 


Some Properties of the Cancellation Multipoles 


Let us consider next some of the restrictions necessarily imposed 
on cancellation-multipole distributions and some of the particular prop- 
erties of those given by equation (36). In the first place, if au(x) 
is any multipole distribution that generates a potential field given by 
equation (8) or (15) , it follows from equation (17) that the value of 
an(x) and its first n derivatives should be everywhere continuous. 
Further, if an(x) is a constant behind some point, say Lq (i.e., for 
oc >x > Lq) , the induced flow field would simulate expanding streamlines 
in the case n = 0 or some form of vorticity in the case n > 0; the for- 
mer case is to be avoided since any simulated body is assumed to have a 
finite area at x = 00 , and the latter case is to be avoided if there are 
no resultant forces normal to the free stream. 


One can show that all the above 
properties are satisfied by an(|), 
the optimum cancellation-multipole 
distribution as defined by equa- 
tion (36). First, notice that 

^n+i)(^) -(jg zero everywhere 

outside the wing-enclosing Mach 
forecone and aftercone, that is for 
-00 < I < -Lo’ and Lq < | < <» (see 
sketch (k)). (Any multipoles in 
these regions cannot combine with 
the wing equivalent multipoles and 
must, therefore. Increase the drag.) 
Hence, one can set 


(n). ^ 

an (-Lq ) 


(n-i) 


(-Lo’) 


= • . . c%(-Lq ) 


r-7x 


tan'/9 

Wing 

plan form 



^ - constant 


= 0 




Sketch (k) 

(n) 

Then the condition of continuity is automatically satisfied for an (|) 
in the entire interval -00 < x < 00 if S^'(|,0) is derived from a wing 
having finite wave drag (in particular, from a wing having no blunt edges 
along which the normal component of the freep-stream Mach number is unity 
or greater). It follows immediately that m < n, is continuous 

since the latter is found by integrating (further smoothing) o^*^^(|). 
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The proof that an(Lo) as given hy equation ( 36 ) is zero requires 
more consideration. One can show, however, that an (Lq) = 0 where 
0 < m < n. First, the equality o 4 "^(Lo) = 0 follows from the fact that 
the wing closes and Sw*(Lo,0) itself is zero. Next consider the defini- 
tion of o4“^(|). Thus 


= ~~ f Sw*(l,e)cos ne de 


cos \x cos ne de / ^ sin cos n)dq 

o ho( 0 ,|) 


- r 


where use has been made of equation (28) for the definition of 
and ho and h^ are defined in sketch (k). Since 


an (I) 


r( 


- 


n-m -1 (n) 

5x) On ^(Si)d|i 


(37) 


one has 


-2Ur 


(-p)“jtr(m - n) 


o2ir 

j ' 


cos ne de 


/ I pB.i(e,li) n-m-i 

d|i / dq(| - li) Xu(^i + ^ sin cos ^)cos 

-Lp* 


Change the li,n coordinates back to the x,y system by means of equa^ 
tion ( 26 ) and this becomes 




-2Vr 


(-p)\r(n-m) 


J' cos ne de J' J' ii- x + cos 


n-m-l 

0) "^y 
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The area S|, shown in sketch (k), becomes independent of 9 when | - Lq 
( being then just the area of the wing itself) , therefore 


4“^Lo) = , .n~ ~r — T f f 

(-13) Jtr(m - n) 


,n-m-i 

X + |3y cos 6 ) cos n 0 d© = 0 


n 0 de = 0 


Hence, for the On(0 defined by equation ( 36 ) 

o 4 ’^^(Lo) = ^^(Lo) = ^^(Lo) = • • • = an(Lo) =0 ( 38 ) 

AIRPLANE SHAPE 


I < 


Lo " 


since, for m < n 


/ 


2Jt 


COS™0 CO 


In the previous section a connection was established between multi- 
pole distributions and their resulting wave drag. Further, this connec- 
tion was direct and relatively simple if the strengths and positions of 
the distributions were given. Unfortunately the connection between the 
multipoles and the shape of the simulated surface is generally not so 
simple. Such a relation does certainly exist, however, and if the 
strengths of the multipoles are known, the relationship is again direct. 
That is, a given distribution of multipoles yields directly, by the for- 
mulas given in the previous section, the induced velocities everywhere 
in the flow field; these, in turn, fix the stream surfaces along any one 
of which (since, of course, the theory neglects viscosity) a physical 
surface can be imagined. 

In general, if 

Fc(x,y,z) = 0 
Fp(x,r,0) = 0 


( 39 ) 
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axe the equations of a stream surface in Cartesian and polar coordinates, 
respectively, then the equations 


BFc SFc ^Fc 

(Uq + Cpx) + Cpy + <Pz = 0 


^ ^Fi, 5Fp 

(Uq + cpx) + 9r + 


1 ^Fp 

x2 he 


= 0 


J 


m 


must hold. 

For example, in studies of thin wings lying in a plane, the partic- 
ular form of equation (39) 


z - h(x,y) = 0 

is assumed and equation (4o) becomes 


-(Uo + cpx) 


Sh 


hh 




» 0 


or, neglecting second-order effects. 


Sx 



which is the familiar boundary condition used in thin-airfoil theory. 
On the other hand, if the equation of the body shape is written in the 
form 


r - R(x,e) = 0 

then equation (40) becomes, for linearized theory. 



= 


R^ he 


(hi) 


If the flow field is radially symmetrical or if the body surface is 
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quasi-cylindrical, equation (4l) reduces to 


Sx 



(42) 


which is the familiar boundary condition used in the study of quasi- 
cylindrical bodies or bodies of revolution. 

In general, a nonlinear partial differential equation of the first 
order such as equation (4l) can be reduced to two simultaneous, ordinary, 
nonlinear differential equations of the first degree (see, e.g., ref. 8). 
Thus equation (4l) becomes 


dx 


1 

UoR^ 


cp0(x,R,0) 


> 


dx 


i cpr(x,R,e) 


( 43 ) 


and if cp^ and are known functions of and 6, these can be solved 

numerically. 

If the strengths of all rectilinear multipoles and source sheets 
are given^ equation (19) or (24) can be used to find cpp and cpg at arbi- 
rary field points. Henc.e^ the first step in finding the body shapes 
reduces to that of integrating such equations. However^ these integra- 
tions are difficult and tedious even when entirely numerical procedures 
are employed and the results still have to be interpreted in terms of 
the body shape according to equations (43). Therefore^ from a practical 
viewpoint^ it is necessary to study certain approximate methods for 
obtaining the velocity field. 


Let attention be concentrated on the disturbances created by a line 
of multipoles. In particular^ consider the fields induced by simple 
polynomial distributions satisfying, in each case, the end conditions given 
by equation (17)* For particular variations set 


ao(50 = Cq(1 

U2(x) = C2(l 
a4(x) = C4(l 


x2)x 

5^2)®x > 

-2x5- 

X ) X J 


( 44 ) 
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where Cq, Cz, an<i C4 are constants determining the amplitudes and 
Ic = x/Lq. Figure l,(a) shows the variations of these coefficients with 
X, and figures 1 (h) throu^ l(f) show how velocities induced by these 
distributions vary with Xo(xq = x - pr) ^d r (r = pr/bo). The results 
have been compared with those for large r given by equation (21) and 
with those for small r given by slender-body theory. Values for the 
latter theory are determined from equations (9) or (I9) by expanding the 
expressions in powers of r and neglecting all but the first terms. 

Thus it can be shown 


ap(x) 

2itr 


^ o = < 


n 


(-2) n!a^(x)cos n6 


4 rtr 


n+i 


n = 0 


n > 0 


( 45 a) 


( -2)'^n!a^(x) sin n6 


r 


r o 


4 irr 


n+i 


n > 0 


( 45 b) 


The significance of figure 1 with regard to practical applications 
is more or less obvious. The first step in its use is to find the effec- 
tive length of the cancellation-multipole distributions. Since the wing 
is given, the streamwise variation of the cancellation multipoles can be 

calculated. Actually this 
variation will extend 
between the apexes of the 
enclosing Mach forecone 
and after cone, a distance 
of Lq + Lp* (see sketch 
(k)). However, depending 
on the wing plan form and 
section, the effective 
lengths of the distribu- 
tions (the interval of 
principal variation) can 
be considerably less as 
illustrated in sketch (Z)). 
Designate this effective 
length as 2Le and the 
distance to the vicinity 
of the body surface as 
rp, and one can define 
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the parameter re thus 


re 





ih6) 


Using figure 1 and the parameter re, one can now estimate the error 
incurred by the use of various approximate methods for calculating the 
body shape. A convenient way to carry out these estimations is to study 
the magnitude of the first crest of the waves shown in figures ,l(b) 
throu^ 1(f), and the distance this crest lies from the foremost Mach 
cone. Graphs showing the variations of these quantities with r are 
given in figures 2 , 3 , and 

By means of the above concepts, let us study briefly four different 
approximate methods that can be used to calculate a body shape. 


Slender -Body Theory 


Slender -body theory is represented in figure 3 by the straight lines 
having the slopes, on the log-log scale, equal to -(n+l) where n is 
the order of the multipole. Since this theory amounts to an expansion 
of the equations for the velocities in powers of f, it obviously repre- 
sents a good approximation when re is sufficiently small. Notice that 
for a given percentage error the limiting value of re for which the 
method applies increases as the order of the multipoles increases. For 
example, when re = 0 - 2 , cPoj. given by slender-body theory is 19 per- 
cent less than that given by exact linearized theory for the case shown, 
whereas cp 2 p is only 3 percent less. Correspondingly, the positions 
of the wave crests follow the path predicted by slender-body theory to 
larger values of r as the order of the multipoles increases. The lat- 
ter trend is illustrated by figure 4. 


If for a particular problem Tq is small enou^ for slender-body 
theory to be considered a good approximation, the equation for the body 
shape, r = R(x,e), corresponding to the combined wing and optimum cancel- 
lation multipoles defined in equation ( 36 ) is determined by the expres- 
sion (using equations (45), ( 36 ), and (l4) together with equation (42)) 


dx 


cpg(x,R,e) 

UqR2 


s’" 

” I f) slnne ^ 2 Jf «x 

VP, 


I 


4n^R 


n +2 


wing 


I dxi(x-Xi) cos nijfSw (x,\|r) 

(47a) 


° -I-o’ 
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dx 


n 


pPj.(x,R,en ^ y 

L Uq -»wing 4r’ 


f i ^ cos n0 

w 

kn E 


.2Jt 


J" dcif J' dxi(x - Xi)^cos n>lfSw”(x,t) 

(47b) 


-Lo’ 


Approximate methods for finding (cpg/UoR^)^j^^g and (<Pr/Uo)ving^ veloc- 
ities induced by the wing, can 
often be used also, hut these 
apply to individual cases and 
cannot he discussed here. 


Source distribution given by 

Off “ CqUoU-x^} X 




Theory for Large Yq 


The asymptotic values for 
magnitude and position of the 
first wave crest obtained by 
placing equations (44) into 
equations (21) are also shown 
in figures, 2, 3^ sind 4. For 
n < 4 it is clear that this 
theory can be used when rg is 
greater than about 2. 

Control-Surface Theory 


The approximations inherent 
in ordinary control-surface the- 
ory can also be estimated by 
inspecting figures 3 4, where 

by control-surface theory one 
means that the exact linearized 
theory is used to evaluate induced 
velocities along a given surface 
and these values are assigned con- 
stant for all r in the vicinity 
of the surface. As shown in 
sketch (m), this amounts to assign- 
ing (Pr and cpg are given by a 
strai^t horizontal line in fig- 
ure 3 and by strai^t lines with 
a unit negative slope in figure 4. 
Obviously, the error in the body 
shape calculated by this theory 
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increases as the amplitudes of the disturbing multipoles increase and as 
the radius of the control surface diminishes. 

One of the simplest applications of control-surface theory arises 
in the study of quasi-cylindrical bodies. In such cases the expression 
for the body surface can be derived immediately from equation (42). Thus, 
if the amplitudes of the cancellation multipoles are small enough and 
the radius of the control surface, is . large enough for control- 
surface theory to be considered a good approximation, the body shape, 
r = R(x,e), corresponding to the combined wing and optimum cancellation 
multipoles is determined by (using equations (I9c) , (l6b), and (36) 
together with equation ( 42) ) 


dx 


Cp (x,R ,e)‘ 


lin ■ 


/ -1 X - i 

0 ° - .x-|3Rc (x-|)cosh(n cosh — 


Z CjjCOs n0 P-^c 

4«2Rc~ J 


-'wing 




*i(x - 




/ S^"U,'lf)cos n\lf d\lf 


(48) 


■where cJn = 1 n = 0 and = 2 for n > 0. 

A study of optimum fuselage shapes using control-surface theory 
has been carried out by Nielsen (ref. 9) for ^ constant-chord sweptback 
wing having a biconvex section and a sonic leading edge. The set of 
interfering singularities used in reference 9 are equivalent (the singu- 
larities ^are limited to the x axis) to the multipoles used herein. 

The fuselage shapes calculated by Nielsen are thus the same - within 
the accuracy of control-surface theory - as those given by equation (48). 


Modified Control-Surface Theory 


A method of modifying control-surface theory to increase its 
accuracy is illustrated in sketch (n). Induced velocities computed by 
this method are based on those calculated along a given control surface 

7 

but extended away from this surface varying in magnitude as (f)- where 
the value of 7^ is fixed by the slope of the curves in figure 3 
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r = re, Tq being defined by equation (46) (see sketch (n)). With this 
modification equation (48) becomes 


dx 


■(Pj.(x,R ,e)' 


Ur 


Ob 

+ 

wing o 



y 

e \ Oncos n6 

T 


4n^Rf 


x-PRq ( X - I ) coshi 


/ 

-Lc 

Source distribution given 

Oq -- c^Uod-x^)x 


(n cosh"^ 

\ PRc/ 




7(x-|)^-p%^ 

by 


J £v"(l,t)cos n\lf d\|r ('49) 


3 . 0 - 


2.5 


2 . 0 - 

/5c^. 
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.15 


.2 


•<c r 
Amplitude of first crest 


.8 


.2 % _.3 

r 

Position of first crest 

Sketch (n) 






1 

1 

1 

1 

Aan'S„ 


1 

1 

1 

1 

1 

1 



4 


which can be solved using 
numerical techniques. 

A further refinement of 
equation (^9) can be obtained 
if the position of the induced 
velocities is also varied 
according to the slope (again 
at r = re) of the curves in 
figure 4. Defining this slope 
as 6n^ see sketch (n)^ and 
xg as 

XB = X - p(l + &n)(R - Rc) 


one can see this refinement 
simply amounts to replacing the 
value of X in the ri^t-hand 
term of equation (49) with xg. 


ILLUSTRATIVE EXAMPLE 
ELLIPTIC WING 


In order that one may be 
able to assess the practical 
significance of the preceding 
sections, the concepts presented 
therein will now be applied to 
the solution of a particular 
problem. For the basic wing 
plan form in this particular 
example an ellipse will be 
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chosen. There are two good reasons for this choice; first, the ellipse 
is effectively unswept and places a severe test on the role of body 
interference in reducing the wave drag at a supersonic speed, and, second, 
for a given volume, the optimum section (i.e., the one yielding minimum 
wave drag) for these wings when considered separately has been discovered 
(see ref. 10) so the reduction in wave drag brou^t about by the body will 
reduce the minimum value possible for such wings when flying alone. The 
drag reductions for the first few cancellation-multipole distributions 
will be calculated and compared with the total drag of the wing alone, the 
wing mounted on an infinite cylinder, and the wing mounted on a basic 
body of revolution. Finally the details of calculating a body shape 
simulated by the wing source sheet, a source line representing a basic 
body of revolution, and the first two optimum cancellation multipole dis- 
tributions will be carried out. 


The Elliptic Wing 


Consider the elliptic lens specified by the equation 


z = ±iri-2e-^' 


b2 


(50) 


where the thickness, span, and chord are shown in sketch (o). The stream- 
wise slope of the upper surface is seen 
to he 


Sx 


u 


= Mx,y) = - ^ (51) 

a 


and the total wing plan-form area S and 
volume V are, respectively. 


S = nab 


V = ^ tab 



(52) 


Wing plan form 


Wave drag . - The wave drag of the 
elliptic wing represented by equation (50) 
can be calculated by means of equation (3O) 
in which, since one wishes now to find the 
wave drag of the wing alone, the an’s are 
set equal to zero. The value of S^'"(x,e) 




2 a- 


Root section 

Sketch (o) 
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follows by placing equation (51) into equation (28) and integrating. Thus 

1/ \ r^V t\ 2tUoCOs / q^sin ^ 

UqS^ (x,e) = 2U0COS [1 J 1(1 +T1 sin ^l)dTl = ^ Uti + ^ J 

^hr 



where, by referring the equation of 
the plan form to the coordi- 

nates (see eq. (26)) and solving 
for the points where the strai^t 
line I = constant intercepts the 


^ wing edges, one finds - see sketch (p) 


hi 

h, 


-b^l sin p ± abJa^cos^p+b^sin^p-l^cos^p 


a^cos^p+b^s in^p 


Sw’(x,e) = - 


l+xtab 


(a^ + b^p^cos^e)' 


Va^ + b^p^cos^e - x2 (53) 


From the relation 

L^(e) = a^ + b^P^cos^e (5^) 

the wave drag can be expressed in the form (integrating once by parts) 

L(e) 


27 t 


o -L(e) -L(e) 


L(e) .,2 

4 tab 


lL^(e)J 


L"( 0 )- 2 |i" 


J'i 


L^(e) 




J'L^ie) - 


l 2 


li- U 


Further integration yields 


q tab 


r — 

J (a^ + 


de 


(a® + b^P^cos^e)' 
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Finally, the wave drag can be expressed in coefficient form, based on 
the total wing area itab, as 




1 + 


2 (“ 

.bp 


1 + 




3/2 


( 55 ) 


Equation (55) represents the lowest value of wave drag possible for a 
wing having an elliptic plan form and fixed volume. This equation was 
first derived by Jones in reference 10. 

The velocities Induced by the wing source sheet in the vicinity of 
the fuselage .- Later, when one wishes to calculate a stream surface in 
the presence of the source sheet that simulates the wing given by equa- 
tion (50) > It is necessary to know the velocities induced by these 
sources at the body surface. Hence, the value of cpr induced by the 
source sheet was calculated at the four points indicated in sketch (q). 
As it turns out, these values are Values of d /U 

so close (see the sketch for a . ^ ^ 

numerical comparison) to those 
obtained by assuming the source 
sheet to be two-dimensional with 
a chordwise intensity identical 
to that along the root section 
of the elliptical sheet (i.-e., 
using the Acker et wave generated 
by the root section) that the 
effect of the wing can be 
to be given everywhere in 
vicinity of the body by the lat- 
ter velocity field if (as will 
be the case in subsequent 
application) the surface of the 
body passes throu^ the region 
shaded in the sketch . That is , 
the effect of the wing in the 
equations for the fuselage shape 
(such as eqs. ( 47 ), ( 48 ), or 
(49)) is assumed to be 
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Sketch (q) 
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1 

rUo 


'Pe 


t / , ■ 

+ ~2 0)sin 0, 

Si 

T-i (x + pr sin 0)cos 0, 
ar 


- for 0 < 0 < rt 
+ for rt < 0 < 2n 

- for 0 < 0 < rt 
+ for rt < 6 < 2rt 


( 56 ) 


The Optimum Cancellation Multipole 


One can now find the strengths of the multipoles along the x axis 
which induce around a cylinder of infinite radius a momentum field iden- 
tical to that created there by the elliptic wing. The negatives of these 
variations are, according to equation (36), the optimum cancellation 
multipoles. Hence, combining equations ( 53 ) (36) 




-On 


2tabUoX N/a^~+~b^P^cos^0~"-~j^ 

Jt(-P)’^'/ (a^ + b^p^cos^O)^ 


cos n0 d0 


( 57 ) 


where On = 1 n = 0 and On = 2 for n > 0 . Particular variations of 

a^^^^(x) are shown in figure 5 - These results are for n = 0 , 2 , and 4 , 

since o^^^(x) for any odd n is zero by symmetry, and apply when the 
wing plan form and free -stream Mach number are related by 


a 

bp " 3rt 


(58) 


which contains the particular case for which the Mach number is n ^2 and 
the aspect ratio is 3. It is apparent that there are at least n + 1 

roots to (Xn^^(x) for -Lq <x < Lq- This follows immediately from 
equation (38) and is true in general. As a result the cxrrves for the 
higher values of n become increasingly wavy and, correspondingly, 
increasingly difficult to evaluate numerically. 

Figure 6 presents the values of a«(x) for the same elliptic-wing 
Mach number relation given by equation (58) . Notice that each of these 
curves has only one root (they necessarily have at least one) in the 
intervail -Lq < x < Lq and is increasingly smooth with increasing n. 

The latter follows from equation ( 37 ) and the fact that the first n 
derivatives of these curves must, in general, be continuous. For example, 
at X = ±Lq the first foiir . derivatives of a4(x) must vanish. 
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Wave drag . - One can now calculate how much the wing-alone drag is 
reduced when combined with each successive optimum cance] lation-multipole 
distribution. If ^ denotes the drag saved by the nth-order cancel- 
lation multipoles, then by equation ( 35 ) 


P' 


Dn 




2 n (n+x) ( n+x) 

J dxi J dxso-n (xx)a-n (x 2 )ln|xi - Xs\ (59) 


-L. 


-L. 


where Lq is the maximum value of L(0) as given by equation (5^) 


= ..2 ^ v2o2 


Lo = a + b"3 


(60) 


The total drag saved by means of the first m multipole distributions, 
would, by equation ( 3 ^) , be 


m 


D - 2Dq + y Dn 

X 


(61) 


,(n+i), 


Using equations (53) and ( 36 ) to define the ''(x) in equation (59) ^ 

reversing the order of integration, and integrating once by parts, one 
finds 


n _li o«/2 oL(ei) pl-iOs) 

— = -1^ / cos n0idei y cos n02d02 / / d^2 

o o _L(03^) _L(eg) 
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32 {tahf 



(a^ 


cos n0x <3-01 
+ b^P^cos^Qi) 



rt /2 


cos nO^dQs 


(62) 


The total drag saved by using all the cancellation multipoles is, by 
definition. 


D 

q 


and since 

00 


n=i 

this is equal to the drag of the wing alone, as it, of course, should be. 

The reduction in wave drag as the wing is combined with the first 
three optimum cancellation multipoles is presented in figure 7. In 
studying figure 7, one sees, as the Mach number approaches one (i.e., 

3 O) , more and more of the original wing wave drag i£ destroyed by a 
line of simple sources alone. Further, the value of re which can be 
written - see equations (46), (54), and figures 5 and 6 - 


-I 


i sin 2nx cos 2nx 


l6(tab)' 
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- = Pre 


(63) 


tends (for a fixed average distance to the body surface re) to zero as 
the Mach number approaches one; and this, in txirn, means that as 3 goes 
to zero the effect of the multipole strengths on the body shape can be 
calculated using slender-body theory. 

When various orders of multipoles are distributed along a line, one 
can show the cross-sect.ional area normal to the free stream of the sim- 
ulated body as given by slender-body theory is a function of the source 
distribution only (see Appendix B) . Coupled with the' discussion in the 
preceding paragraph, this can be used to demonstrate that, for Mach 
numbers close to one, the "supersonic area rule" proposed in reference 11 
and discussed in reference 12 gives a good approximation for the wave 
drag of an elliptic wing and body combination which is symmetrical with 
respect to the plane of the wing. 

The induced velocity field . - A method for calculating the velocity 
field induced by the multipoles when On(x) is given numerically is pre- 
sented in Appendix C. By means of this method, velocities induced by 
the do and cxz multipole distributions shown in figure 6 have been 
calculated for r equal to 0.li;8 and the results are shown in figure 8. 
Since the distributions; in figure 6 were for the particular case 
a/b3 = k/ 3 n, it is evident from equation (63) that the values in figure 8 
apply to the case Tq/Io equal to O.I6I; that is, when the body radius 
is about 16 percent of the wing semispan. 

For comparative purposes, the values given by slender-body theory 
are also shown in figure 8. The degree of agreement between the two 
cu3Tves is consistent with the results shown in figures 3 


Interpretation of Drag Reductions 


Comparison with wing moiinted on a 
circular cylinder .- With regard to fig- 
ure 7 oae should be careful to notice 
that the drag of the wing alone has been 
used for the reference drag. The drag 
reductions shown, therefore, represent 
gains brought about by interfering with 
the velocity field induced by a planar 
source sheet, or, in tenns of a combina- 
tion with an upstream cylindrical stream 
surface, gains made by modifying a body, 
shown in sketch ( s) , which bulges behind 


Mach waves from wing 
roof section 



\ \ 


Sketch (s) 
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the wing leading-edge Mach wave in accordance with the velocities induced 
there by the source sheet. Obviously, from this viewpoint, a considerable 
reduction in drag can be brou^t about merely by eliminating the bulge, 
thereby making the body a circular cylinder throughout. Mathematically, 
such a procedure amounts to using a certain set® of cancellation multi- 
poles along the x axis behind the point -Lq, and, if the drag of this 
resulting combination were used as a reference, the gains shown in fig- 
ure 7 would be diminished. . 

An approximate way to estimate the drag of a wing mounted on a cir- 
cular cylinder is illustrated in sketch (t) and consists merely of sub- 
tracting from the wing source sheet those sources blanketed by the body. 



Sketch (t) 


Using the subscripts 1, 2, and 3 to designate the wave drags of the indi- 
vidual wings as indicated in the sketch, Jones (in an unpublished com- 
munication) has shown that if wing 1 is an elliptic wing with a biconvex 
section and wing 2 lies entirely within the plan form of wing 1, then 



m 


where V 2 and Vi are the volumes of wings 1 and 2, respectively. 


0 

The exact evaluation of multipole distributions necessary to simu- 
late a circular cylinder for the entire body length has been studied in 
reference I 3 . 
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For a supersonic Mach number, Dg is closely approximated by the 
wave drag of a rectangular wing having the same section and aspect ratio. 
If Ag, ^Ha^ and Tg are, respectively, the aspect ratio, plan-form area, 
and thickness ratio of the rectangular Ving, its drag can be expressed in 
the form 


^ ^ l^RaUg (65) 

qT2 


where 


Ng = 



3Ag 


2 sin~^Ag 

3 3Ag 



16 

3 


PAg<l 

(66) 


PAg>l 


Further, if A^ is the aspect ratio of the elliptic wing, one can show - 
see sketch (t) and equation (52) - 


The drag of the elliptic wing follows from equation (55) and can be 
written 


^ = i^abNi (68) 

qxg 


where 



( 69 ) 


Finally^ therefore^ equation (64) can he put in the form 
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^ = fi . ^ K 

L 3it V'b/J b Ni 


(70) 


and the ratio N 2 /N 1 is a function 
By means of equation (TO)^ the 


of the parameters R/b and ^A^ only, 
dashed curves shown in sketch (u) — 



Oo 




representing approximately the wave drag of a wing mounted centrally on 
a circular cylinder - were calculated. Though considerable drag reduc- 
tion is indicated by adding just those multipoles necessary to make the 
body cylindrical, it is apparent the. total wave drag can be reduced fur- 
ther, for the range of parameters shown, by using only the first two 
optimum cancellation-multipole distributions, aQ(x) and a 2 (x), given by 
equation ( 36) . 

Comparison with wing mounted on a basic body of revolution . - 
Sketch (u) shows the effect on the wave drag of adding the optimum can- 
cellation multipoles either to the wing alone or to the combination of 
an infinite circular cylinder and a centrally mounted wing. Estimates 
of their effect when added to a wing mounted on a basic body of revolu- 
tion can also be carried out. In order to present these estimates, how- 
ever, the results of the following two theorems due to R. T. Jones are 
needed. 
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1. Designate the closed body of revolution which, by a slender- 
body theory, has a minimum drag for a fixed volume and length as a Sears- 
Haack body. Then the total wave drag of a Sears-Haack body and any other 
body of revolution or any centrally mounted thin xkng which lie entirely 
within the Sears-Haack body's enclosing Mach forecone and aftercone is 
given by the equation 



+ Dg 


(71) 


wave drag of Sears-Haack body alone 
wave drag of other body or (exposed) wing alone 
Vgg volume of Sears-Haack body 
V 2 volume of other body or (exposed) wing 

2. Designate the closed body which, by slender-body theory, has a 
minimum drag for a fixed base diameter and length as a Karman ogive. 

Then the total wave- drag of a Karman ogive and any other slender body of 
revolution or any centrally mounted thin wing which lie entirely within 
the ogive’s enclosing Mach forecone and aftercOne is given by the equation 


where : 

dsh 

D2 


D — Dj^ + D2 


(72) 


where : 

Dg wave drag of K^m^ ogive alone 

D 2 wave drag of other body or (exposed) wing alone 

In order that the theoretical results could be tested by wind-tunnel 
experiments, a basic body of revolution having a finite base area was 
chosen. Such a body can be simulated by a combination of the source dis- 
tributions which produce,® separately, 

B ■■■ ■ - 

The source distributions simulate the Sears-Haack body and the 
Karm^ ogive only when slender-body theory is used to calculate the body 
shapes. If linearized theory is used instead, the body shape will, of 
course, be somewhat different. However, the subsequent results and con- 
clusions are by no means limited by the assumptions pertaining to slender- 
body theory. The latter theory is used only to obtain an estimate of 
the body volumes or to study cases for which it gives results that are 
not significantly different from those given by linearized theory. 
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the Sears-Haack body and the Kannan ogive. Thus, if 21 is to be the 
body length, the line of sources 

^o(x) 2 x\J 




2 

itl3 


Vk - hY 


SH 


,2 2 

1 - X 




o Jll- \ V 

simulates (by slender-body theory) 
a body of revolution (see sketch (v)) 
having a total volume V equal to 
VgH + cross-sectional area given 

by 


Sketch (v) 

[.J7'- 

and a base area S(Z) equal to Vj^/Z 


Vk 


Z^(^ + sin“^y 


n 8V, 




3/2 


3«Z 


-Z < X < Z 


(7l^) 


The wave drag of a wing mounted on this basic, unmodified body will 
now be calculated. Just as was the case in studying the wing attached to 
an infinite cylinder, the assumption is made that the wave drag of this 
combination is the same as the wave drag on the configuration simulated 
by superimposing the singularity distributions which create separately 
the exposed wing panels and the body of revolution. With this assumption, 
the wave drag can be written explicity in terms of the wing and body 
geometry by applying equations (71) and (72). Hence, 


D = % + D| 




+ D. 


(75) 


where D3 is the drag of the exposed wing panels alone ^ given by equa- 
tion (70) and shown for various values of R/b in sketch (u) , and V3 
is their volume (see sketch (t)) 


IT 4 . 'uf ^ 8 R 

V3 = tatj{ 5 ^ 


( 76 ) 


Since D^; and Dgg, the wave drags of a Karm^ ogive and a Sears-Haack 
body flying alone, are well known to be 




« Z‘ 


8q VsH 

»SH = - — 

the wave drag coefficient of the unmodified combination, based on the 
complete wing area nab, can be expressed as 


(77) 


qnab 


= Cp = 


Vjf + 8v 


SH 


l6tv 


SH 


n^Z'^ab 


n^Z^ 


8 £ - 

3 b/ na^p 


t^Ni 


32 r\ R £2 
3« V b Ni. 


(78) 
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where Ni and Ng are defined in terms of Mach ntunher and wing-body 
geometry in equations ( 66) and ( 69 ) . An example of the variation of Cp 
with Mach number for the particular combination shown in sketch ( x) 

(R/b was set equal to 0 .I 8 I) is given by the dashed line in sketch (w). 
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It is now possible to find how much the drag of this unmodified com- 
bination can be reduced by means of the optimum cancellation-multipole 
distributions used to derive the results shown in figure 7 . Again apply- 
ing equations ( 7 I) and ( 72) , one can show 
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D = Dk + Dgn^ + + ^4 ( 79 ) 


:i_, 


- where V4 is the total voliane and D4 
the total wave drag pertaining to the 
wing (now the complete wing including 
the portion blanketed by the body) and 
the multipoles. However, within the 

Maximum thickness of wing along accuracy of the approximation - being. 


center tine , t , - 0. 234 

Total volume of body - 44.60 
VsH- 3 1.72 




- 12.88 


Sketch (x) 

a specific value of | 3 Ai 
can readily show 


in fact, exact within slender -body 
theory, see Appendix B - the volume 
added by the wing is subtracted from 
the basic body by the optimum cancel- 
lation - source distribution so that 
V4 is zero. Further, if N4 is the 
value of D/Dw read from figure 7 for 
and a specific number of multipole types, one 


D4 = 


t^q 
3 a 

where is defined in equation (69) 

combination can be reduced to either 


NiN. 


1«4 


(80) 


Hence, the drag of the unmodified 


qitab 


= Cp = 


n^abZ 


|vk 


+ 8 


V, 


SH 




jta^ 


NiN, 


liM4 


(81) 


if the same total volume is maintained (maintained, as is obvious from 
an inspection of the equation, by increasing the value of the Sears- 
Haack portion of the basic body an amount equal to the volume of the 
exposed wing) or to 


CD 




jta^ 


NtN, 


11H4 


(82) 


if the volume of the fuselage is reduced by an amount equal to the wing 
volume . 

The results expressed by equations (80) and (8I) , when applied to 
the first two optimum cancellation-multipole distributions, are shown 
for the geometrical parameters presented in sketch (x) by the solid 
curves in sketch (w) . The value of R/b used for the solid curves was 
0.161 instead of the O.I8I value used to calculate the dashed curve. 

The smaller value was used since the modified body is drawn in along the 
sides by the cancellation multipoles (see fig. 9 ), decreasing the average 
body radius in the wing region from about 1.00 to about 0.89. One must 
be careful to notice that the solid curves represent minimum (relative 
to the special method being discussed) values which can be obtained by 
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a specific design at a specific Mach number and do not represent the 
variation of wave drag with Mach number for any given combination. 


The Body Shape - First Calculation 


The final step in studying the effect of the optimum cancellation 
multipoles , defined in equation ( 36 ) , is to find the distorted body shape 
which they produce in combination with the wing and a basic body. The 
decision was made to calculate a body shape which would be optimum at a 
Mach number equal to ^/^. The details of the wing and body geometry are 
given in sketch (x) and the basic body parameters Vk and Vg were inter- 
preted in terms of source strength by equation ( 73 )- 


It was apparent from the results of figure 8 that, for the values of 
rg and b given by sketch ( x) , the velocity field induced by the first 
two optimum multipole distributions can be calculated with good accuracy 
using slender-body theory. Combining the values of Tj. and opg so cal- 
culated with those induced by the wing, given by equation ( 56 ), and those 
induced by the basic body, using also slender-body theory to interpret 
equation ( 73) , one can find the body shape by solving the two simultaneous 
nonlinear differential equations presented as equations ( 43) . These were 
solved numerically by the method outlined in Appendix D and the results 
were, unfortunately, unrealistic. Sketch (y) shows an example of a 


streamline close to the 0=0 plane 
and the crossing of such streamlines 
obviously invalidates the solution. 


Body Shape - Second Calculation 


Ihe failure observed in the 
first calculation has a simple enou^ 
interpretation. For the chosen wing 
the basic body was too small in diam- 
eter at the wing-body juncture to 


Wing 

source 

sheet- 



Muifipoies 


'Streamline^ 0-^0 

Sketch (y) 


permit the use of the first two cancellation multipoles in their entirety. 


Several avenues of approach are yet available. One could, for 
example, maintain the same wing and basic body but reduce the Mach number, 
one could start with a larger basic body, or one could lower the thickness 
ratio or aspect ratio of the wing, thus diminishing the strength of the 
cancellation multipoles. All of these, however, are modifications of the 
basic conditions or basic restraints and as soon as such restrictions are 
abandoned it must be remarked that no matter how low the' wave drag of a 
set of nonlifting, volume -enclosing siarfaces has been made, another 
arrangement of the same volume within a finite space will give a lower 
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vfti IIP j unless^ of course^ tlie wave dz*ag of the first arrangement is 
already zero. Therefore, instead of modifying any of the initial 
restraints, consider the following alternative: 

How much can the drag he reduced by using only a portion 
of the first two cancellation multipoles so that a real 
body would still be simulated? 


In order to answer the above question, examine briefly the first 
calculation. Notice, from figure 8, that the cause of the body collapse 
is attributable to the large values of q)j. and cpg induced by the second- 
order cancellation multipoles, clz{x) . Hence, let ao(x) be maintained at 
its full value and reductions permitted only in the magnitude of a 2 (x). 

To cairy out such a procedure efficiently, one must be able to determine 
the effect of a given variation of tt 2 (x) on the body shape. Fortunately, 
Graham, in reference l4, has developed a method by means of which the 
relation between ct, 2 (x) and body shape can be quickly estimated. Graham 
has shown, if 


M -- /.4/ 



1. a rectilinear distribution of 

second-order multipoles of 
strength a^ix) is placed 
along the x axis (see 
sketch (z)) in a supersonic 
stream (M = \l2 ) , 

2. slender-body theory is used 

to evaluate cpj. and cpg, 

3. R = Rq is the radius of a 

circularly cylindrical tube 
for -Lq > X > - 00 , 

-Lo 


then 


R = Rc 


1 + 2f’^(x)cos 20 + ^/l + 4F^(x)cos 20 + 4(F*(x)]‘ 




(83) 


is the continuation of the stream tube for x > -Lq. 


Since the initial strengths of the cancellation multipoles are 
negative, F*(x) is negative and the critical value of R occi^s along 
the plane 0 = 0 or n. The v^iatlon of R/Rq with, is given 
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in sketch (aa) from which one can see that the maximum value of If*| 
must be less than 0.5 if the simulated stream tube behind the plane 
X = -Lq is to be real. 

The problem can now be continued, using Graham’s result as a guide, 
by assuming the critical body radius in the more complicated source and 
multipole arrangement is principally determined by F*(0) - the parameter 
governing the body indentation at the center of the cutout and at the 
wing-body juncture. In the first 


place, since there is no interference 
between different orders of multipoles 
it is necessary to consider only the 
drag produced by the second-order 
multipole. Appendix E presents a 
method for finding the optimum dis- 
tribution of the second-order can- 
cellation multipoles for a given 
wing and a fixed value of F*(0) . 

The resulting wave drag is given in 
equation (Ell) . At a Mach number 
equal to \/2 and for the basic wing 
and body parameters presented in 
sketch (x), F^*(0) (defined by eq. E8) 
equals - 2.90 and the reduction in 
Da/qS, the amount of drag caused by 
the wing second-order multipoles 
alone - see equation ( 62 ) and fig- 
ure 7 ~ is shown in sketch (ab) for 
a range of F*(0). Variations on 
the strength of various combinations 
of second-order multipoles are shown 
in sketch (ac) on page il8. It is 
important to notice that for a given 
percentage reduction in the maximum 
strength of the multipoles th£ result 
ing percentage reduction in Dg is 
much larger. 

The strengths of a 2 (x) shown 
in sketch (ac) must now be combined 
with the zero-order multipoles and 
wing source sheet, and the combined 
velocity field used to calculate the 
shape of the new body. Using again 
slender-body theory to evaluate the 
velocity field Induced by the multi- 
poles and the numerical methods given 
in Appendix D to compute the stream- 
lines, one finds, by restricting the 



Sketch (aa) 



Sketch (ab) 
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Sketch ( ac) 

distribution and strength of the second-order cancellation multipoles to 
their optimum values corresponding to the restraint F*(0) = -0.6, 
that a rea l as well as reasonable body shape results. The details of 

i°Sketch (aa) gives 0.5 as the maximum permissible value of -F*(x). 
However, that value is based on a distribution of Oe(x) alone, ^d in 
our more general case the added velocity field caused by the presence 
of the other singularities permits the larger value. 
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this shape are presented in figiare 9 and their general interpretation 
is discussed in the next section. Finally, using the value F (0) = -0.6, 
the drag curves shown in sketch (w) were reinterpreted, and the results - 
which represent an estimate of the amount the wave drag of an elliptic 
wing mounted on a basic body of revolution can be reduced by realistic body 
distortions - are shown in sketch (ad) . 
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Discussion of Results 


9x 


It is possible to gain some insist into the reasons for the various 
body distortions shown in figure 9 by inspecting, in another light, the 
body shape first calculated. Consider first the elliptic wing at the top 
of sketch (ae). The air over the forward part of this wing, when it is 
alone in a supersonic stream, is compressed (mathematically, the sign of 
is negative), the compression being greatest near the leading edge. 

On the other hand, the air over the after portion 

I Uo of the wing is undergoing an expansion, the magni- 

tude of which is greatest near the trailing edge. 
Consider now, in combination with this wing, a body 
which is to have a shape providing favorable inter- 
ference. It is apparent that the body should cast 
expansion waves over the forward portion of the wing, 
destroying the compression there, and absorb the 
expansion waves coming from the wing after portion. 

Or in another li^t, the positive pressure on the 
forward region of the wing (one can use the equation 
Op = (p-Po)/q. = "2q)x/Uo for the pressure coeffi- 
cient) should be reduced as far as possible by a 
wave shed from the body and having large negative 
pressures where it comes in contact with the wing 
forward region. 

Since waves in a supersonic flow field are 
fundamentally associated with the slope of the 
disturbing surface, the aforementioned favorable 
interference fields would be created by a body 
having, longitudinally along its surface, slopes 
such as those shown in the lower part of sketch (ae). 
This is exactly what the solution obtained from the 
calculation of the first body shape tried to estab- 
lish since the fuselage near the plane of the wing 
(the portion most strongly affecting and being 
affected by the pressures on the wing) and ahead 
of the wing chordwise center line was distorted 
in a manner that caused an expansion across the 
wing entire forward portion. The difficulty arose 
because the fuselage was not wide enou^ to provide 
the longitudinal extent of favorable slopes neces- 
sary to create the positive pressure called for by 
the wing forward compression region, and the body 
streamline near the wing root, following a path 
such as that shown by the line .in sketch (ae), 
crossed the body center line before it reached the 
wing chordwise center. 



expansion 


Sketch (ae) 
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Consider now the second body calculated in the previous section. 

In this case an additional restraint was imposed which, effectively, 
fixed the maximum body indentation. Subject to such a condition, an 
optimum interference field was discovered. If the resulting fuselage 
shape is inspected near the plane of the wing, surface slopes are found 
similar to those shown in sketch (af). The following discussion is 
intended to show that, from a physical viewpoint, this arrangement is 
reasonable. 


Most of the wing pressure drag occurs on the 
Hence, for a fixed maximum fuselage indentation, 
an' over-all point of view, to create initially a 
increases the pressure drag on the forward por- 
tion of the wing tip but provides a succeeding 
extent of fuselage having slopes that generate 
a strong expansion wave over the forward por- 
tion of the wing inboard section. Similarly, 
the final portion of the body is forced to have 
a region of u nf avorable interferences ■where the 
expansion waves from the wing tips combine with 
body expansion waves to increase the local drag 
(i.e., increase the local suction pressure) in 
order that the over-all interference effects 
are as beneficial, under the given restraints, 
as possible. This arrangement (i.e., unfavor- 
able interference near the wing tip and favor- 
able interference near the ■wing root) is given 
further support by the attenuation property 
inherent in three-dimensional waves. Thus the 
pressures induced by the body ©n the ■wing tips 
are not as strong, for a given generating sur- 
face slope, as those induced on the inner por- 
tion of the wing, simply because the tips are 
farther from the disturbing surface. 


wing inboard portions, 
it is beneficial, from 
compression wave, which 


Body 

sfream/ine 



— Muffipo/es 


Region of 
compression 


expansion 


Although these considerations are somewhat Sketch (af) 

oversimplified (the shape of the upper part of the body has been com- 
pletely ignored in estimating the effect of the waves) , the longitudinal 
variation of surface slopes near the plane of the wing and the resulting 
body streamlines there are, from a physical point of view, reasonable. 


In order to support the above conclusions, the source and multipole 
distributions simulating the final modified body shown in figure 9 were 
used to calculate (see Appendix C) u/Uq in the plane of the wing near 
the root section. The values of u/Uq induced by the wing sources along 
these sections were assumed to be the same as those induced by a ■two- 
dimensional biconvex section having the same local chord; that is, tip 
effects were neglected. These values for body and wing were added and 
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the resulting pressure distributions, shown in sketch (ag) (Cp = -2 u/Uq) , 
were obtained. The results sne similar to the estimates presented in 
sketch ( af ) . The large drag saving near the root section is illustrated 




in sketch (ag) by the graph showing the low values of section drag coef- 
ficient along the inner portion of the wing. 

Another important characteristic of wing-body combinations designed 
to have low wave drag is also illustrated in sketch (ag) . As shown 
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in the graph of u/Uq^ over the surface of a two-dimensional biconvex 
section the air is everywhere accelerating in the streamwise direction. 

In studies concerning the effects of viscosity on the fluid flow and, 
in particular, studies concerning the boundary layer, this positive fluid 
acceleration is referred to as a favorable pressure gradient. If the flow 
is laminar in the vicinity of the leading edge of a smooth wing and the 
pressure gradient is everywhere favorable, the flow tends to remain lam- 
inar and \mseparated over most of the wing chord. Notice that the modi- 
fied wing-body combination has a line of zero pressure gradient extending 
along a Mach line downstream from a point near the body and wing leading- 
edge j\oncture. Immediately behind this line the pressure gradient is 
unfavorable which gives rise to the possibility of flow separation or, 
at least, transition from laminar to turbulent flow there. 


Comparison With Experiment 


The modified wing-body combination shown in figure 9 was tested in 
the Ames 2- by 2-foot transonic wind tunnel . The Reynolds number of the 
test, based on the mean aerodynamic chord, was approximately 1.5x10®. 

This combination had an exposed wing volume of 3-^^ cubic inches and a 
body volume equal to l)-4.60 cubic inches, for a total volume of 48.04 
cubic inches. As a control, an unmodified combination composed of the 
same elliptic wing mounted on a body of revolution (the area distribution 
of which was determined from equation (7^) with Z = 10.5# = 12.88, 

and Vg£ = 29.02 cubic inches) was tested. The exposed wing area in the 
unmodified combination was 3.32 cubic inches and the body volume was 
41.90 cubic inches, for a total volume of 45.22 cubic inches. Thus, the 
unmodified combination had the same body length as the modified one but 
less volmne. 

. The wave drag at M = 1.4l of the combination shown in figure 9 has 
edready been calculated and presented in sketch (ad) by the curve pertain- 
ing to real body shapes. By use, in equation ( 78 )# of the values of Vk 

VsH mentioned above and a value of O.I 76 for R/b, the wave drag 
for the unmodified body was calculated throu^out a supersonic Mach num- 
ber range. The theoretical results obtained for body configurations are 
shown by the dashed curves in sketch ( ah) . 

The wind-tunnel results for the total drag on both configurations 
are shown in figure 10 for 0.6 <M < 1.^1-. Notice that three groups of 
data are shown. The lower one represents the unmodified body alone, the 
middle one represents the modified and unmodified combinations with no 
fixed transition, and the upper one represents both combinations with 
transition fixed along the leading edge. The models tested with natural 
transition did not show the predicted drag reduction. As was pointed out 
in the discussion of sketch (ag), however, the adverse pressure gradients 
on the modified model could be inducing transition in the vicinity of • 
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the 45 ° line behind the juncture of the wing leading edge with the fuse- 
lage. This, in turn, would cause the wing of the modified model to have 
a larger area covered with a turbulent boundary layer and, hence, cause 
the drag of the model to Increase. In order to separate the potential 
and viscous effects, the transition-fixed tests were made. If the exper- 
imental wave drag is taken to be the difference between the drag at a 
supersonic Mach number and the drag at M = 0.6, the resulting values of 
experimental wave drag are as shown in sketch (ah). 
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Sketch (ah) shows that the experimental reduction in wave drag 
brought about by the modification agrees with that predicted by theory. 
Both theory and experiment show a reduction of about 0.0015 in the drag 
coefficient at the design Mach number (l.4l), and the experiment further 
shows an average reduction of 0.0020 over the Mach number range 
1.2 < M < l.U. A further study of figure 10 shows that the difference 
between the experimental and theoretical wave drags shown in sketch (ah) 
for the wing-body combinations is nearly the same as the difference 
between experiment and theory for the body alone. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Feb. l6, 1955 
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APPENDIX A 

DERIVATION OF THE OPERATIONAL FORM OF THE WAVE EQUATION 


For convenience, take the normalized form of the wave equation in 
Cartesian coordinates, thus 


Q 

Sz^ 

and define the Laplace transform of cp(x,y,z) by 


(Al) 


9(s,y,z 



®^cp(x,y,z)dx 


(A2) 


Now if X = f(y,z) is the equation of the foremost Mach cone or Mach cone 
envelope and f(y,z) >0, it is apparent 



f cp(x,y,z)e"®^dx 

oy ^ 



Sy ® W/x=f 


(A3) 


since (cp)x=f 
tion (A3) we see 



is not necessarily zero. From equa- 



o 


Sy 


Sz 




(Al^) 


(A5) 


Further, integrating by parts gives 


r g-sx ^ _ 

J ox^ 



(A6) 
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Hence, 



o 


_ 

S z^J 


dx 


2 - 

■ ^z2 - 

-sf/^ + ^ ^ ^ 

® \Sx By By Bz 


(A?) 


The last term on the ri^t is the directional derivative of the per- 
turbation potential along the surface x = f(y,z). This is, of course, 
along the so-called conormal. Since cp is a constant on the forward 
envelope, its gradient along the envelope is zero and 
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APPENDIX B 

ON THE VOLUME OF BODIES CALCULATED USING SLENDER -BODY THEORY 


The following proof shows that in a rectilinear distribution of 
singularities, only the sources contribute to the total cross-sectional 
area of the simulated body and, hence, to its volime. 

According to slender-body theory, the velocities induced in the 
field by distributions of multipoles along the x axis can be written 


ap(x) 

2itr 


= S 


n 


(-2) n!an(x)cos n© 


4itr 


n+i 


n = 0 


n > 0 


(Bla) 



(“2)^n!an(x)sin n© 


4itr 


n+i 


(Bib) 


Further we have derived - see equation (4l) - neglecting only second- 
order effects, the equation representing the boundary condition for the 
body, thus 


Sr 
° Sx 



*^© Sr\ 

^®/r=R(x,©) 


Combine equations (Bl) and (B2) 


2«UoE |S = a„(x) || 

1 


d©. 


(B2) 


multiply by 


and integrate 
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SJtUo ^J' ~ <ie = 2jtao(x) (-2)“n!an(x)y' 


2JT/R COS n0 - sin n0 

de 

n+i 


de 


or 


00 

2«Uo|^S(x) = 2rtao(x) + ^(-2)“(n-l)!a^(x)y'^'' de (B 3 ) 


Since the integrand in equation (B3) is a periodic function in 6 , ve 
have 


M 

ax Uo 


(B4) 


which shows the simulated-body normal cross-sectional area to be dependent 
only on the source strength. Further, the total volume is given by 



x)s'(x)dx + (Z* + Z)S(-Z) 


and when S(l') =S(-Z) =0, there results 


V 




xeio(x)dx 


(B5) 


-I 
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APPENDIX C 


ON TTni: CALCULATION OF VELOCITIES INDUCED BY ARBITRARY 
SOURCE DISTRIBUTIONS 


The potential and velocity fields represented by equations (l8) and 
flQ) are difficult to evaluate analytically even if An(l) is a simple 
function. However, the calculations can be reduced to a relatively simple 
process. First, let equation (l8) be expressed in terms of the dimension- 
less variables x,l, and f where 


X = x/Lq ' 
I = I/Lq 
f = pr/Lo 


(Cl) 


Then 


x-r 


Ao(i)dg 

5 - i)^ - 


(P2(x,r,e) = " 


2jt 


I 

-1 


X-r 


cos 26 r Ag(|)^2(x - l) ~ r jd^ 
r2 J{x - l)^ - 


(C2a) 


(C2b) 



and so forth. Consider next 
the variation of An(i) shown 
in sketch (ai) and represented 
by the equation 


And) 


Sketch (ai) 

m —2 

( 190n - 360m) + ( 390n - 760m) | + (200n - UOOm) | , 


n 


-1 < f < -0.9 
-0.9 < I 


(C3) 

As seen in the sketch. An vanishes at | = -1, is a 

and -0.9 (assuming the values m at 1 - -0.95 n at | - *9)^ 
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the strai^t line, An(l) = n, for I greater than -O.9. The velocities 
induced by a multipole distribution given by equation (C 3 ) can be cal- 
culated in a strai^tforward manner in the two regions - 1 + f <x < -O.9 + 
and -0.9+ r<x. For example, if 


L ( X + 0 . 9)\/(x + 1.0)^-r^+ f ^In 


X + 1.0 + '/ (x + 1.0) 


2 -2 
- r 




jtfMQ_(x) . 

— — ^ = ^-(x + 0.9)/(x+1.0)^ -f^+ (x+1.0)7(x + 0,9)^-f^ + 


x+ 1.0+V (x+ 1.0)^ - 

r'^ln , -0.9 + r<x 

x + 0 . 9 +nT(x + 0 . 9)2 - f 2 

(C 4 ) 


jtfNop(x) 


t 1 ■ x + 1.0 + /(x + 1 . 0 )^ - f' 

(S +0.95)>/x+ 1 . 0 )^ - f ^ - f^ln 1 


100 


- 


-1 + f <x <-0.9 + f 

(x + 0.95) 7(x+ 1.0)^ -f^ - (x + 1.05) /(x + 0.9)^ -r^ - 

-0.9+ r < X 


r^in 


X + 


1.0+n/(x+ 1.0)^ -r^ 


X + 0.9 +/(x + 0.9) 


2 n2 


(C 5 ) 


then cpQ- can be written 


cpof = mMoj;(x) + nNoji(x) 


(C6) 


Now, if one is given a distribution of sources that is composed of, 
or is adequately approximated by, a series of 20 equally spaced parabolic 
arcs, equation (C6) can be used for each individual arc and the results 
superimposed for the complete solution. To this effect, define mj^ and 
ni in terms of Ao(x) by 
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"1 ■ ■ 0 ■ ■ 0 


(C7) 


so they represent the magnitudes shown in sketch (aj). Then, if [x] 
denotes the greatest interger contained in x (e.g. , [ 6 - 3 ^] equals 6 ), 



[ lo(x-f)+ll] 


90J 


i=i 


(c8) 


Values of Mq- and Nq- are ta^bulated in table I for f equal to 

0 . 074 , 0.l48, 0.222, and O. 296 . The asymptotic magnitudes of these func- 
tions are given by slender-body theory. Hence, one can easily show for 
large x 


M 


Of 


0 



(C9) 


Notice that both functions have essentially reached their asymptotic 
values for large x by the time x=-l+f+0.5. By applying simple 
tabulative procedures to equation (C 8 ) - for example, listing m^ and n^ 
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in reverse order and accumulating multiplications of adjacent terms - 
the value of cpoj; for any Aq(x) representable by equation (C7) is 
readily calculated. 

The velocities induced by higher order multipoles can be calculated 
in a similar fashion. Because of the asymptotic behavior of the M^s and 
N^s^ however^ one is led into the numerically inefficient process of 
obtaining small numbers from differences of large numbers. For the veloc- 
ities cpp and cpg^ the following is a method for circumventing this dif- 
ficulty. 

It follovs from equations (9) or (I 9 ), that for small r, cpnj and cpng 

can be expressed in terms of the multipole strengths an(x) - as defined 
by equation (7) - by the equations 

n+i/ \ 

^ + . . . + C^r^ + . . . (clOa) 


n+: 


«Pn 


6 


r sin ne 


= Dq + D^r + 


+ 


(ClOb) 


where for v < n, n > 1 


■ D (V 


Cv = < 




r 1 + 


(x), V even (Clla) 


V odd 


nr 

(-1) ^pVgn-vhF 






kn 


r(i + I 


(v), . 

an (x), V even 


, V odd 


(Cllb) 


Consider now the velocities cpg- and 
strength defined by equation (C3); thus 


induced by the multipole 
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<P2- 

— = mM2c(x) + nN2=(x) 

cos 29 ^ ^ 


(C12a) 


<Pa 


9 


r sin 29 


= wMsq(x) + dN2q(x) 


(C12b) 


where values of the M*s and N*s are listed in tables II and III. 
asymptotic values, as given by equations (CIO) and (Cll) , are 


19 + 20x 

= li" 

^ 150nf® 


Their 


No- = 


150P + 275x + 126 


150jtf ® 


,, 19 + 2QX 

M 2 _ " _ o 

e I50rtr® 






l^Ox^ + 2T^x + 126 X 

150rtf® 


(CI3) 


and these are also given in the tables. 

As the tables show, equations (CI3) are sufficiently accurate 
approximations to M and N for practical calculating purposes when 
x>-l+r+ 0.5. Hence, referring to sketch (ak) one can see that the 


Effect of multipoles in 
this interval on velocities 
at (x/) given by equa~ 
tion (Ci3h 


/ (x,r) 



/ Line of 
y muitipoies 


L 




f 


-/ 


x-r-5 


x-r 


Sketch (ak) 

velocities at the point x,f induced by the multipoles in the interval 
-1<| <x - f - O.5 can be calculated using equations (C13)» In terms 
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- ( 2) - 

of the distribution for - which is equal to see equa- 

tion (16) - this means the multipole distribution shown in the upper 
part of sketch (az) can be calculated by means of the asymptotic formulas 
and the result added to that obtained for the distribution shown in the 
lower part of sketch (al) by use of equations (C12b) and tables II and 
III in a manner identical to the one represented by equations (C6), (C7), 
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where P^a 2 ^^(x) 
by the equation 


Asix) and is shown in sketch (aZ) and defined 

[10(5 - r)] 


x±=0.h + x- r- 


10 


(C15) 


the symbol [ 10(x - r)] meaning, as before, the highest integer value 
contained in 10(x - f ) . A similar result can be derived for (pgj; and 
one has finally for -1 < [10(x - f) + 11] < 5 


q>2e(x) 
f sin 29 


[ lo(x-f) + ii] 




+ 


n^Ns 




(Cl6a) 


Cp2i:(x) 

cos 29 


[ lo(x-f )+ll] 



1 



+ 



and for [ io(x - r) + ii] > 6 



(Cl6b) 


920(x) 

f sin 29 


[io(x-?)+ii] . 



[ io(x-f)+ii]-5 


+ H(x,f) 

(ClTa) 


92r(x) 
cos 29 


[ io(x-f)+ii] 

I , 

[ lo(x-f)+ll] -5 





Aa(x - f - xj) 
2nf 


(C17b) 


The streamwise gradients of induced velocities can also be defined 
in terms of M*s and as were the velocities cpj and i cp^. Thus 
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«Pox = + nNox 


cos 29 


- wMs- + nNsx 


Values of Mqj, Nqj, Mgj, and Ng- for f equal to 0.IU8, 0.222, and 
0.296 are given in tatles IV and V. 
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APPENDIX D 

NUMERICAL METHOD USED TO CALCULATE BODY SHAPE 


The method used to calculate the body shape was a standard step- 
by-step solution to the two simultaneous total differential equations 
(eq. (i^3) in the text) 




UoR" 


^ i q>r(x,R,e) 


dx U| 


(Dl) 


The essentials of the process are recognized from the following computing 
sheet heading set up for initial values of 6 and R equal to 30 and 
0.ll»8, respectively, where x = x/Lq and R = PRA'o 


Q 


© 


© 

© 

© 

n 

X 

^ XO.05 + 

n-i 

e 

^ xo.05 + 

n 

n-1 

R 

fl 

VUo 8 e/ 
for 

Qand0 

f 

© 

f 1 

VUo 8 fy 
for 

(3) and® 

0 

-0.852 

0.524 

0.148 

- - - 

- - - 

- - - 

1 

1 

bo 

0 

ro 

- - - 

- - - 

- - “ 

“ “ - 

- - - 

2 

-r ?52 

- - - 

- - - 

- - " 

- - - 

- - - 


® 

® 

@ 

© 

1 M 

© 

0 

XO.O5+ 

n-i 

n -1 

61 

® X0.05+ 

n-i 

© 

n -1 

Ri 

{1 M 
VUo 8 e/ 

for 

®and( 9 ) 

©. 

© 

fl ScpV 
VUo 

for 

®and® 

/^de'N 

Vdx/ 

^ av 

© +© 

fd^ 

Wav 

© +@ 
n -1 n 

2 

2 

- - “ 

- - - 

- - - 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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APPENDIX E 


o 

OPTIMUM variation OF ac(x) FOR A FIXED VALUE OF J ct 2 (x)dx 


•“Lr 


Given 


Dp = - 


q r \ (. 3 ) ,(3), ,]r(3), , (3), ^ 

iwu”^ J J (Xl)J &2 (X2) +02 (X2) 

“Lo “Lq 


] 


In 


Xi - X 2 


dxidx2 


(El) 


(3) 


where th« variation of a2 (x) is fixed, pose the restraint 


/ o 

ae(x)dx = F (O) ^ R^* = constant 


(E2) 


-Lr 


and ask for the function 02(x) which minimizes D2 for a given value 
of the constant. 

If f(x) = f(-x), then 


.1 


o o 


J' J' f(xi)f(x2)ln|xi - X2|dxidx2 = 2 y' J' f(xx)f(x2)ln| x^^ - X2^|dxidx; 


■i -1 


•1 -1 


and since 02(-Lo) = oi^\-Lo) = 02 ^\-Lo) = 0 


y 02(x)dx = - x^o^ \x)dx 


-Lo 


Therefore, the standard variational problem 
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Dg + A 


/ 


° a 2 (x)dx 

~ 


= 0 


(E 3 ) 


reduces to 


J J ai®^ (xi) + 02^^(xi)j jai^^(x2) + a!^\^z) 


2 itUo' 


-Ln “Lq 


In 


Lo^ I 


Xi‘"-Xo‘^, , X 


y'xx®o 4 ^\xi)dx 


= 0 


-Lr 


and this becomes 


r\ 

r y1 

tt2 (Xx) 

dxi ' 

q r° 

J ® 

Uo 

«Uo J 

0 

1 

_ -J 


-Lo 


( 3 ), , ( 3 ), ; 

E 2 (X2) + tt 2 (X2) 


In 


Y 2 Y 2 

Xj X 2 


dX 2 + 


Axi^ 

6 


= 0 


Integrating three times by parts, using the relations 


6 a 2 (-Lo) = 5 ai^^-Lo) = 5 ai^^-Lo) = 0 

( 2) 

8a2(0) = 8a2 (O) = 0 , by symmetry 


and 


lim P a2®^(x2) +a2®^(x2) 

xx » o ^ L 


In 


2 2 

Xx - X2 


dxs 


= lim 2xx / 

Xx-^O J 




dX 2 = 0 


Xx - X2 
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yields 


-Lo -Lo 


(3). . (3) . 

as (xa) + as (xs) 


In 


Y 2 Y 2 

Xx -X 2 


dx2+ UoAoXi' 


= 0 


where 


■V _ 


(e4) 


By the fundamental lemma of the calculus of variations 


Sx£" 


/ 

-Lf 


( 2) / \ ( 2) / \ 

aa (xa) +os (xs) 


In 


- X2^ 


Lo‘ 


dX2 +U 0 A 0 X 1 " 


= 0^ 0 > Xx > -Lc 


One can also show 


axx" 




(3) / X (3) / X 

aa (X 2 ) +os (xs) 


In 


2 2 
-Xa 


T 2 


dx2 -UqAoXx' 


= 0, Lq > Xx > 0 


Hence , 


r SS-iMln 
J Uo 


2 2 
Xi -X2 


T 2 




rl 2 -x 3 r ^2 (^ 2 ) 

v7 uq 


In 


V 2 „ 2 

Xi -X 2 ' 




dxs 


(E5) 


Integrating by parts and changing the notation so that 


H 2 - ^ 9 

Lo 


y ^*(^ 2 ) = 


^ , 8 ( 12 ) = 


Bl “ T 2 ^ 

■*->o 


( 2 ), X 

aa ( Xa) 
aa^^(x2) 


one has 


f(B2)d'n2 -X 

/ ~r ;r- = Aq + AaBx “ ^oBi 

J Til - B2 ° 


3/2 g(T|2)<lTl2 

J Tlx - Tia 


(E6) 
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Equation (E6) is the familiar sin^lar integral equation known^ in 
aerodynamic applications^ as the airfoil equation. Its inversion is 
discussed^ for example^ in reference 3* one solves equation (e 6) and 
applies the conditions 


then 


4^^(0) = ai^^Lo) = 0 
15rt ^ 


and 


(XaCx) ^ ag(x) XqLq^x 
Uq Uq 60 jt^ 



(E7) 


Wow set 



ag(x) 

Uo 


dx = F/(0) 


(E8) 


so that 


Fw*(0) 


is a known constant. 


Then 


CX.2(x)dx 


= F* (0) I Rd*= -Fw*(0) J Re** 


Ur 


2AqLo^ 

225rt^ 


-Lr 


(E9) 


Using the above expressions, one can show 


°^g(x) 

Uo 


a2(x) 


Ur 


19«Rc^x 
32L 2 


* 

F (0) + 




(ElO) 
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The wave drag can be calculated by combining equations (E5) and (El) 




aa (xi) aa (xi) 
- + 


Uo 


Ur 


dx- 


Integrate three times by parts and there results, finally 


B _ 675 . 2 3 ^^ 


F*(0) + F/(0) 


(Ell) 
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Figiare 1.- Continued. 
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Figure 1.- Continued. 
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Figure 1.- Continued. 
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(e) Velocities, cp 4 ^. 
Figure 1.- Continued, 
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Figure 1.- Continued. 
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Figure 1.- Continued. 
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(f) Concluded. 


Figure 1.- Concluded. 
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(b) Radial velocities induced by second-order multipoles. 
Figin:e 2.- Continued, 
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(c) Tangential velocities induced by second-order multipoles. 

Figure 2.- Continued. 
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( S’) <Por 

Figure 3 . - Attenuation of first crest of waves created by multipole dis 

tributions shown in figure 1 . 
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Figure 3.- Contiriued. 
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Figiire 3-“ Concluded. 
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(c) Tangential velocities induced by second-order multipoles. 

Figure !<•.- Continued. 
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(d) Radial velocities induced by fourth-order multipoles 



(e) Tangential velocities induced by foiirth-order multipoles 

Figure 4 .- Concluded. 
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(a) n = 0 

Figure 5*“ Variation of nth derivative of nth-order cancellation multi- 
poles for elliptic wing. 
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Figiire 5-" Continued. 
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(a) n = 2 

Figure 6.- Variation of nth-order cancellation multipoles for elliptic 

wing. 
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(b) n = If- 

Figure 6.- Concluded. 
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(a) Radial velocities induced "by sources. 

Figure 8.- Velocities induced by the elliptic-wing cancellation multipoles at the control surface 

where PR/Lq = 0.lil8. 
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(b) Radial velocities induced by second-order multipoles. 
Figure 8.- Continued. 
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Figure 10.- Experimental values for drag coefficient. 
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